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Abstract 

Spontaneous chiral symmetry breaking in QCD with massless quarks at infinite volume 
can be seen in a finite box by studying, for instance, the dependence of the chiral con- 
densate from the volume and the quark mass. We perform a feasibility study of this 
program by computing the quark condensate on the lattice in the quenched approxi- 
mation of QCD at small quark masses. We carry out simulations in various topological 
sectors of the theory at several volumes, quark masses and lattice spacings by employing 
fermions with an exact chiral symmetry, and we focus on observables which are infrared 
stable and free from mass-dependent ultraviolet divergences. The numerical calculation 
is carried out with an exact variance-reduction technique, which is designed to be par- 
ticularly efficient when spontaneous symmetry breaking is at work in generating a few 
very small low-lying eigenvalues of the Dirac operator. The finite-size scaling behaviour 
of the condensate in the topological sectors considered agrees, within our statistical ac- 
curacy, with the expectations of the chiral effective theory. Close to the chiral limit we 
observe a detailed agreement with the first Leutwyler-Smilga sum rule. By comparing 
the mass, the volume and the topology dependence of our results with the predictions 
of the chiral effective theory, we extract the corresponding low-energy constant. 



1 Introduction 



Spontaneous symmetry breaking plays a central role in our understanding of the strong 
interactions. In QCD with (a small number of) Nf light flavors the standard expec- 
tation is that the SU{Nf)L x SU{Nf)ji chiral symmetry group breaks spontaneously 
to SU{Nf)L^pi with the formation of a non-zero quark condensate. There are many 
evidences that support this picture [1-14], but a conclusive study with a reliable deter- 
mination of the condensate is still missing^. 

The lattice formulation of gauge theories is at present the only approach where 
QCD can be defined non-perturbatively. Computations in this regime are performed by 
numerical simulations, and therefore limited to be at finite volume and lattice spacing. 
The spontaneous formation of a quark condensate S at infinite volume can be detected 
in a finite box by studying the properties of the long-range correlations left in the 
system, a technique widely used in statistical mechanics [16]. In QCD the simplest 
observable to consider is maybe the chiral condensate generated by adding a small 
quark mass m to the action. This explicit breaking of chiral symmetry gives also rise to 
additive ultraviolet divergences which make the definition of the condensate ambiguous. 
The crucial observation here is that the ultraviolet divergences are independent on 
the volume V and, in a regularization that preserve chiral symmetry, their magnitude 
is suppressed by a factor 1/V with respect to the contributions from the long-range 
correlations [17]. At any given lattice spacing the value of S can then be extracted with 
arbitrary precision from a finite-size scaling study of the condensate if the volumes are 
large enough. 

In a finite box it is also interesting to consider QCD in a fixed topological sector, 
i.e. with the functional integral restricted to gauge-field configurations with a given 
topological charge. Recent progress in the understanding of topology on the lattice [18- 
20] allows to show that the chiral condensate has the same ultraviolet divergences in 
the full theory and at fixed topology. It can be made ultraviolet finite with the very 
same renormalization constants and subtraction coefficients. On the other hand the 
chiral effective theory makes definite predictions, with no extra free parameters, for the 
topology dependence of the infrared contribution to condensate [17,21,22], predictions 
which can be directly tested against lattice QCD results. 

Despite of many efforts dedicated to implement this program in the lattice com- 
munity [7,10,11], a reliable determination of the QCD chiral condensate is still missing. 
The main reasons being the limitations in our ability of going beyond the quenched 
approximation in the numerical simulations of QCD, and the large numerical instabil- 
ities in the commonly used estimators of the chiral condensate at finite volume [11]. 
The first problem is being solved over the last few years thanks to the development 
of new algorithms for simulating dynamical fermions [23,24]. The implementation of 
some of these new ideas to simulate dynamical Ginsparg-Wilson fermions is progressing 

^For a recent review of the experimental implications of spontaneous symmetry breaking in QCD 
see Ref. [15] and references therein. 
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quite fast. First results were already presented [25-28]. The second problem arises 
in computations of correlation functions at finite volume in the mass range rriY^V^l. 
The long-range correlations left in the system generate small eigenvalues of the Dirac 
operator of order 1/V [17]. The latter give rise to extremely large fluctuations in the 
commonly used statistical estimators of correlation functions, and make it difficult to 
control the statistical errors [11,29]. 

In this paper we address the second problem. We introduce numerical estima- 
tors of the condensate which are infrared stable and free from mass-dependent additive 
ultraviolet divergences. This progress builds on the recent understanding of the renor- 
malization properties of the spectral density of the Dirac operator [30]. The form of 
these estimators represents and explicit example of how the informations from the chiral 
effective theory can be used not only for the interpretation/extrapolation of the data, 
but also to design more efficient numerical algorithms. 

We apply these ideas to the computation of the chiral condensate in the quenched 
approximation of QCD with Neuberger fermions. For cubic lattices with linear exten- 
sions L^l.5 fm, we observe a finite-size scaling behaviour of the condensate compat- 
ible with a theory which exhibits spontaneous symmetry breaking at asymptotically 
large volumes. The topology dependence of the data supports the predictions of the 
(quenched) effective theory. Close to the chiral limit we observe a detailed agreement 
with the first Leutwyler-Smilga sum rule [17]. Eventually we extract the value of the 
low-energy constant S by matching the chiral effective theory formulas with the lattice 
results. 

The paper is organized as follows. In Sec. [2] we discuss the ultraviolet properties of 
the chiral condensate with Neuberger fermions; in Sec. [3] we summarize the predictions 
of the effective theory for its volume, mass and topology dependence; in Sec. [4] we 
give details of the numerical computation, and in Sec. [5] we compare our results with 
the predictions of the effective theory. The appendices \M and [B] are devoted to more 
technical details. 

2 The quark condensate with Ginsparg— Wilson fermions 

We consider QCD with Nf > 2 degenerate flavours regularized on an Euclidean lat- 
tice of spacing a with the standard plaquette gauge action and with the Neuberger- 
Dirac operator [31]. The latter (see Appendix A for unexplained notations) satisfies the 
Ginsparg- Wilson relation [32] 

75D + Z)75 = aDj^D , (2.1) 

which guarantees an exact chiral symmetry 

§ip = (1 — aD)ip , 511) = ipj^i , (2-2) 

of the massless fermion action at finite lattice spacing [33] . The Jacobian of the trans- 
formation is non-trivial, and the chiral anomaly is recovered d la Fujikwa [33-35], with 
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the topological charge density operator defined as [31,36,37] 



a\{x) = --iTYi^D{x,x)\. (2.3) 
In a given gauge background the topological charge operator 

Q = a''Y.q{x) (2.4) 

X 

is related to the index of the lattice Dirac operator Q = n+ — n_ where (n_) is 
the number of zero modes of D with positive (negative) chirality. All the cumulants of 
the topological charge Q are ultraviolet finite and thus the charge distribution [18-20]. 
Correlation functions of renormalized local operators inserted at a physical distance are 
therefore ultraviolet finite also in the theory restricted to a fixed topological sector. The 
action of the theory is invariant under the non-singlet chiral transformations analogous 
to the one in Eq. (|2.2p . They ensure that in the chiral limit the quark condensate 
renormalizes only multiplicatively 

(V^V) = lim Zsii^i^) , (2.5) 

where ijj = {V'l, . . . , V'TV/}) is defined analogously and the rotated field V'i is given in 
the Appendix A. Zs is the logarithmic-divergent renormalization constant of the scalar 
density fixed with a given renormalization prescription. For our finite-size scaling study 
we are interested in breaking the chiral symmetry explicitly by adding a mass term to 
the action, which can be introduced by defining the massive Dirac operator as 



D. 



m 



am\ „ 

l- — ]D + m. (2.6) 



2 J 



The quark condensate develops additional ultraviolet divergences which can be param- 
eterized as _ _ 

_ _|_ ^^777,3 _|_ {finite terms} , (2-7) 

with the asymptotic behaviour of bi and 62 being l/o? and ln(a) respectively. In finite 
volume it is useful to define the chiral condensate in fixed topological sectors v = \Q\. 
It can be parameterized as 

where 

Xi/ = ZsXu = bim + b2m^ + {finite terms} . (2-9) 

The 1/m infrared divergence in Eq. ()2.8p is the trivial topology-dependent contribution 
due to the presence of the zero modes of D in the quark propagator. The renormalization 
constant Zs, and the additive subtractions, which are needed to remove the additive 
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ultraviolet divergences, can be chosen to be independent on the topology, see Appendix 
|BJ As a consequence the continuum limit of the combinations 

Xu,-Xu, (2.10) 

is unambiguously defined at finite quark mass provided the volume is large enough. 

The spectral density p(A) of the massive Dirac operator Dm has a well defined 
thermodynamic limit. It can be defined as [38] 

/'W = ^E(^(^-i^^i)) (2-11) 

k 

where |Afcp are the non topological eigenvalues of D'^D restricted to one of the chiral 
sectors. Recently it has been shown that its renormalized counterpart 

p{X) = Zsp{Zs\) (2.12) 

has a universal continuum limit [30]. The same applies to the spectral density Pui^) 
defined in the theory at fixed topology. At this point the renormalized spectral density 
p{X) can, in principle, be computed with lattice simulations. However it is more realistic 
from a numerical point of view to consider integrals of /5(A) with a given probe function. 
There is clearly a lot of choice. In the rest of our paper we will focus on observables 
"condensate-like" of the form 

r.(A^in, A„,ax) = 2m I \ f ^ ^ p,{\) dX , (2.13) 

-'An.in [1 - (aZs m/2)2]A2 + 

with A = X/Zs and rh = m/Zs, which have a well defined continuum limit if Amin 
and Amax are kept fixed when a — > 0. For Amax Aqcd these integrals are expected 
to be dominated from the most infrared part of the Dirac spectrum (see below) where 
also discretization effects are expected to be quite small. A more detail study of other 
possible choices is left for future studies. At fixed lattice spacing the quantity 

oo) = ^ lim T;^ (Amin, Amax) (2-14) 



has the same ultraviolet divergences of the chiral condensate in Eq. (j2.9p . They arise 
from the contributions to the integral at the upper end of the integration domain. Since 
they are topology independent, the observables 

T!/i(Amml, oo) - f^2(Amin2,Oo) (2.15) 

are unambiguously defined. The chiral condensate is recovered by taking Amin = 
in Eq. ()2.14p . The coefficient of the leading order in m is the first Leutwyler-Smilga 
sum rule [17] 



1 A 

V dm 



Xu 



r '-^™" M^)^^. (2.10) 

Jo X^ 



_ 2^ 1 - {aZs X/2y 

m=0 y Jo A2 

and differences of sum rules analogous to Eq. (j2.10p are free from ultraviolet ambiguities. 
All these considerations remain valid in the so-called quenched approximation, where 
the fermion determinant is dropped in the effective gluon action. 
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3 The quark condensate in the effective chiral theory 



In presence of spontaneous symmetry breaking, the QCD correlation functions at small 
masses and momenta are expected to match those of a chiral effective theory which 
contains only the pseudo-Goldstone bosons as dynamical degrees of freedom [1,2]. At 
asymptotically large volumes L ^> 1/ Aqcd and for small quark masses the finite volume 
behaviour of QCD correlation functions can then be predicted [39,40]. In the mass range 
rhTiV^ 1 the partition function of the effective theory reads 

Z = [ dC/o exp [^ ReTr t/o] , (3.1) 

JsU(Nf) 

where C/q are the zero-momentum modes in the chiral theory, dUo is the corresponding 
Haar measure, fi = rhT,V, and S is the infinite volume chiral condensate. The quark 
condensate (with opposite sign) is given by [17] 

S(/z) = -^^InZ = -^/i- <( ■' +••• , (3.2) 




Nf djj. 2Nf 



and it vanishes in the chiral limit since spontaneous symmetries breaking can only occur 
in infinite volume. For fi <^ 1 the linear proportionality of the condensate to the quark 
mass and the volume is a distinctive sign of a theory which undergoes spontaneous sym- 
metry breaking in infinite volume [16,17,38]. It is due to the long-range correlations left 
in the system at finite volume. The quark condensate plays the role of the magnetization 
of a ferromagnet with the quark mass being the analogous of the external magnetic field. 
A tiny variation of the quark mass requires a variation of order (mV) in the condensate. 
This behaviour is not spoiled in QCD at finite lattice spacing if fermions preserve an 
exact chiral symmetry. The mass-dependent ultraviolet divergences in Eq. (j2.7p affect 
only the sub-leading corrections suppressed by a factor 1/V. 

The partition function of the theory at fixed topological charge z/ at asymptotically 
large volumes and in the same mass range reads 

Z^= [ (detC/)^exp[/xReTrC/] , (3.3) 

Ju{Nf) 

and the quark condensate is given by [17] 

MM) = 5f + X.W = ^|;l..Z. = ^ + 5(^M+-- A'/>2 (3.4) 



The infrared divergence matches exactly the analogous one on the QCD side in Eq. 
For ^ <C 1 the linear proportionality of the condensate to the quark mass and the volume 
signals in this case the presence of a non-zero condensate in the thermodynamic limit of 
the full theory. The low-energy constant S can be determined with arbitrary precision, 
if the quark mass and the volume V are properly tuned, by matching the combinations 

XuAf^) - XuiifJ-) (3-5) 
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with the QCD counterparts in Eq. (j2.10p . The next-to-leading (NLO) corrections to the 
scahng behaviour can be computed in so-called e-expansion [40]. They are suppressed 
as {AttFL)~'^ , with F being the Goldstone boson decay constant in the chiral limit, and 
they can be included by replacing S with the volume dependent parameter TiesiV), 
whose expression can be found in Ref. [40]. This implies that the ratios 



~ 7 \~t — I AT \t — rTr\ ' '^s 7^ , (•3-Dj 

^=0 {^■i-^i){vi + Nf){v2 + Nf) 



are non-trivial parameter free predictions of the effective theory up to NLO, which is the 
higher order at which sub-leading corrections are presently known. The functional form 
of the condensate is known for arbitrary values of /i < 1 [17,41]. The condensate 

is also known in a partially quenched chiral effective theory with Nf flavors and an extra 
valence one of mass rhy. From its discontinuity across the imaginary axis 



Disc 



,SPQ(/i, m,) = hm SPQ(/i, i\ + e)- SpQ(^, i\ + e) = 27rp,(A) , (3.7) 



it is possible to extract the NLO effective theory prediction p(A) of the spectral density 
in presence of Nf flavors [21,22]. This is the functional form that is expected to match 
/5(A) in Eq. (j2.12p for A ^ Aqcd ^ 1/a. In the following for brevity we report the 
full expressions for the chiral condensate and the spectral density for the quenched case 
only, which is the one needed in rest of the paper. 

In the quenched approximation of QCD, an effective low-energy chiral theory is 
formally obtained if an additional expansion in 1/Nc, where Nc is the number of colors, 
is carried out together with the usual one in quark masses and momenta [42,43]. Here 
we adopt the pragmatic assumption that the quenched chiral theory describes the low- 
energy regime of quenched QCD in certain ranges of kinematic scales at fixed Nc- 
Correlation functions can be parameterized in terms of effective coupling constants, the 
latter being defined as the couplings of the effective theory. The quenched condensate 
at fixed topology at asymptotically large volumes (LO) is given by [21,22] 

= M(/,(/.)K,(^) +W(/.)K,_i(/.)) + -, (3.8) 

where 1^, are modified Bessel functions. At the NLO in the e-expansion the con- 
densate is [21,44] 

with /ieff = mTiefiV. The infrared sickness of quenched chiral perturbation theory, as 
well as of quenched QCD, manifests itself in the fact that diverges logarithmically 
with L. In the following we define operatively at a given volume as 



rhV 



= ^i^{V) u>Q. (3.10) 

m=0 
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Figure 1: The microscopic spectral density and PuiOK"^ foi' ^ = 0^5 are shown on the 
left and on the right plots respectively. 



If we consider two box sizes Vi^2 = L\ 2: the volume dependence of at NLO in 
e-expansion of the quenched chiral theory is [44,45] 

2 1 

^2 



where /3i = 0.140461 is a shape coefficient [46], tjiq and a are the singlet mass and 
the additional low-energy constant of the singlet kinetic term in the chiral Lagrangian 
respectively (see Ref. [44] for unexplained notation). From the discussion above it follows 
that the ratios in Eq. (j3.6p . for i/j > and with Nf = 0, are non-trivial parameter- 
free predictions also in the quenched effective theory up to the sub-leading corrections 
of higher order presently known. The spectral density of the Dirac operator in the 
quenched approximation can be extracted from Eqs. (j3.7p . (13. Sp and (j3.9p . They lead 
to the microscopic spectral density [21,44] 

PuiO = ^ [UCf - J.+i(C) J.-i(C)] , (3.12) 

which is related to as 

PuiO = ^uiC/^V) . (3.13) 

with ( = XTiV. For A <^ Aqcd the integral Ti,(Amm) Amax) in Eq. ()2.13p can then be 
matched with 

Cmin Cmax\ _ ^_ p^^^^ rf^ ^ (3.14) 
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The functions Pu{C) PviOK"^ ^'^^ shown in Fig.[TJ They highhghts the fact that large 
parts of the integrals in Eq. (j3.14p come from the region with C~20, where a handful 
of low-lying eigenvalues of the Dirac operator provides the bulk of the contribution. 
Similar conclusions apply in the range < 1. 



4 Numerical computation 

The primary observables in our numerical computations are the integrals ^^{Xjam^ oo) 
defined in Eq. ()2.14p . The lower bound Amin is chosen to ensure the reliability of the 
Monte Carlo estimate (see below) and at the same time to minimize the difference with 
Xy. The low energy constant Sgfr is extracted from a matching of the combinations in 
Eq. ()2.15p with the corresponding formulas in the effective theory. Our final goal, the 
computation of {xu^ — Xv2 ) ■> is then achieved by adding the small contribution from the 
infrared tail analytically using the functional form in Eq. (|3.14p . 



4.1 Numerical estimator 

The underlying chiral symmetry guarantees that for every background gauge configu- 
ration the quark propagator Sm{x,y) (see Appendix [A]) satisfies 

Y,tT[s^{x,x)] = - + 2^tr[p,S„(x,x)P,] (4.1) 

X ^ X 

with v being the absolute value of the topological charge of the configuration and Pc 
the chiral projector into the chiral sector without zero modej^. In this form, and once 
the topological charge is known, the calculation of the quark condensate requires the 
computation of the propagator in the sector without zero modes only [7]. 
On reasonably large lattices the computation of all eigenvalues of the Dirac operator 
or of the propagator from every source point x is numerically unfeasible. It is quite 
standard, however, to compute the quark propagator from a few fixed source points or 
to extract a few low-lying eigenvectors. Let us assume that once has been determined, 
a number n of approximate low-lying eigenvalues of PcD^DPc and the corresponding 
eigenvectors is computed by minimizing the Ritz functional starting from random vec- 
tors generated with a Gaussian action, i.e. invariant under space-time translations [47]. 
The minimization is carried out until the approximated eigenvalues have a relative error 
(jjk and satisfy 

PcDWPcUk = \Xkfuk + rk, {k = l,...,n) ; (4.2) 
{ui,uk) = dik ; (4.3) 
{ui,rk) = Oy{l,k), \\rk\\ < uoklXkf . (4.4) 



In finite volume the probability of having a configuration with zero modes in both chiralities is zero. 
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The propagator in the sector without zero modes can then be spht in a light and a heavy 
contribution as follows: 



^ PcUk{x)ul{y)Pc 
k=i (1 -a^"^V4)|Afc|2 + m2 



n 



where Uk is defined from as in Eq. (jA.Op . It is easy to prove that, once averaged over 
the gauge configurations, the spin-color trace of each contribution on the right-hand side 
of Eq. (j4.5p is translational invariant even if the Uk are only approximate eigenvectors, 
i.e. Wfc 7^ [47]. The condensate, after the trivial contribution from the zero modes is 
subtracted, can then be decomposed as 

X- = -^E(t4'^-(^'^)l) -^=xl + X^ (4.6) 
where the heavy and the light contributions can be computed as 

= ^-"^I^^IV^ \ (47) 

X^l = 2 (tr[p,5^0,0)pJ) . (4.8) 



It must be stressed that Eq. (j4.6p is exact independently on the values of ujk and the 
number of the extracted eigenvectors n. By contrast, the statistical variance of the esti- 
mator changes with oj^ and n. The local fluctuations of the approximated eigenvectors 
Uk{x) are enhanced on the r.h.s. of Eq. (j4.5p by the smallness of the denominator. Un- 
der the working assumption that they are responsible for large variations in the trace 
of the local propagator, the variance of the estimator in Eq. (j4.6p is greatly reduced 
with respect to the one of tr [Pc'S'(0, 0)Pc]- Deep in the chiral regime the probability of 
having a configuration with eigenvalues lying in the infrared tail of the spectral density 
/5(A) is small but not negligible. The tail is poorly sampled by the Monte Carlo while at 
the same time the smallness of the eigenvalues generates large fluctuations in the light 
part of the condensate xL- overcame this problem we replace with 

in Eq. (|4.6p . where 6 is the usual step function. The value of Xmin is chosen in such a 
way that only 5 eigenvalues in all the Monte Carlo history of a given data set are lower. 
The parameter SgfT is then extracted by matching these observables with the analogous 
ones in the effective theory. Eventually the chiral condensate Xv is computed by adding 
the small contribution from the tail analytically using the formula in Eq. (j3.14p . When 
the latter is substantial the determination of Egg from differences of f,y(Amin) oo) is still 
correct, but the value of Xv is heavily affected by the functional form used and therefore 
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lat 


f3 L/a L (fm) iVcfg 




am 


cl 


5.8458 12 1.49 fm 672 


119, 205, 155, 
104, 51, 29 


0.001, 0.003, 0.008, 
0.012 ,0.0016 


c2 


5.8458 16 1.98 fm 488 


49, 69, 82, 
72, 50, 54 


0.000316, 0.000949, 0.00253, 
0.00380, 0.00506 


c3 


6.0 16 1.49 fm 418 


74, 137, 101, 
62, 27, 12 


0.000612, 0.00184, 0.00490, 
0.00735, 0.00980 



Table 1: Parameters of the simulations: /3 = 6/gf^ is the bare gauge coupling, L is the 
linear extent of the each lattice, A'^cfg is the total number of configurations generated, 
and am are the bare quark masses considered. N^^^ refers to the subset of configurations 
with fixed topological charge for z/ = 0-5. 



it less interesting to us. In the following we only consider data sets where the integral of 
the tail is at most 25 per cent of Xu, and in most of the cases it is less than 10 per cent. 
A positive side effect of this requirement is that the computation of the eigenvalues does 
not have to be very precise. The tail contribution in Xu can, of course, be reduced by 
increasing the statistics of the data set, i.e. sampling part of the tail with confidence. 

This procedure complements for the case of the chiral condensate the low-mode 
averaging (LMA) technique proposed in Ref. [48]. The latter has been successfully 
applied already to meson two-point functions [48-50], baryon two-point functions [51] 
and more recently to three-point functions for the extraction of low-energy constants of 
the AS = 1 chiral effective Hamiltonian [52]. 

4.2 Numerical experience 

We have generated ensembles of gauge configurations by standard Monte Carlo tech- 
niques with the Wilson gluon action and periodic boundary conditions. The fermions 
are discretized with the Neuberger-Dirac operator as defined in Eq. (|A.3p with s = 0.4. 
A summary of the parameters of our runs are reported in Tab. [TJ We have simulated 
three lattices with two different volumes (runs cl and c2) and two lattice spacings (runs 
cl and c3). The linear extent always satisfies L^l.5 fm, a size that we expect to be large 
enough for a finite-size scaling study. This is suggested by the results in Ref. [12], where 
for L^l.5 fm it was found a detailed agreement of the predictions of random matrix the- 
ory for the low-lying eigenvalues of the Dirac operator with quenched QCD results. The 
values of the two lattice spacings are chosen to guarantee the locality of the Neuberger 
operator [53], and to be in a range where discretization effects were found to be small in 
several observables [13,54]. The quark masses for the lattice cl are fixed to be roughly 
in the interval 0.07 ^ /-i ^ 1.2. For the lattice c2 the masses are such that the values of 
(mV) match those of cl. The masses for the lattice c3 are chosen so that the dimension- 
less quantity (mV/ ZsTq) is constant, where vq = 0.5 fm is a widely used reference scale 
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Figure 2: Monte Carlo history for the run cl, am = 0.008 and for the topofogical sectors 
ly = (left) and i/ = 1 (right). The first plot from the bottom shows the Monte Carlo 
history of the absolute value of the smallest eigenvalue of PcD^DPc, the second one 
represents the condensate (divided by the quark mass) computed without LMA (black- 
thin line) and the sum {tI + Xu) / i^''^) (magenta-thick line). The third and the fourth 
show the Monte Carlo histories of the heavy contribution {a^x'^) / am and the light one 
{a^rl) / am, respectively. 



in quenched QCD computations [55, 56] and Zs is the renormalization constant of the 
scalar density in the RGI scheme which we have taken from Ref. [13]. The calculation 
of the topological charge, of the low-lying eigenvalues, and of the quark propagator is 
performed following Ref. [47]. For each run the low- mode averaging is implemented as 
described in the previous subsection with n = 20 and uo^ = 0.05. A posteriori we have 
verified that the highest eigenvalue extracted satisfies ((|A2o|)Sefrl^) > 20 for all lattices. 

In Fig. [2] we show a typical Monte Carlo history for the smallest eigenvalue of 
the Neuberger operator, for the chiral condensate with and without LMA and for the 
heavy and light contributions separately. It is obtained from the run cl for i/ = 0, 1 
and am = 0.008. A first observation is that the heavy part of the condensate is very 
stable in all topological sectors. Moreover the quantity tr PcS^{Q,G)Pc /{am) is essen- 
tially independent on the quark mass configuration by configuration. For all lattices, its 
largest relative deviation that we have observed among different masses is roughly 10~^. 
We interpret this as a consequence of the fact that |A2o| ^ m configuration by config- 
uration. We thus expect that a further stabilization of the heavy contribution would 
not reduce the variance of the condensate significantly. The Monte Carlo history of 
does not show large spikes for all masses, and a statistical analysis is applicable. Some 
of the lighter masses at lower topologies, however, have been discarded (see Tab. [2]) to 
satisfy the upper limit on the tail contribution discussed in the previous subsection. As 
expected the light part fluctuates much more than the heavy one. Large contributions 
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appear in coincidence with the lower values of |Af | consistently with the expectations. 
The Monte Carlo history of the local estimator tr [PcS{0, 0)Pc] has fluctuations which 
are much larger than those with LMA. Moreover for the lower topologies extreme sta- 
tistical fluctuations are observed for most of the masses considered, which invalidate the 
statistical analysis of the sample. In the following the physics analysis is carried out 
only on the data with LMA. 
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Figure 3: The quantities a'^xt/'>T^ (bottom) and a^xl (top) as a function of u, for the 
second heaviest mass of the run cl (black points) and run c2 (red points). 

In Fig. [3] we show our results for the heavy and the light contributions to the chiral 
condensate as a function of for the second heaviest mass of the runs cl and c2. The 
heavy part divided by the quark mass is weakly dependent on the topological charge 
and the volume. The results for a^xL &t the corresponding masses are compatible. 
This indicates that the light contribution is to a good approximation a function of the 
variable (mV). In addition, the 1/z^ behavior predicted close to the chiral limit by the 
first Leutwyler-Smilga sum rule is reproduced qualitatively. These observations, which 
are valid for all masses simulated, point to the fact that the splitting between heavy 
and light contribution is such that the bulk of the heavy part behaves essentially like an 
ultraviolet divergence (i.e. volume and topology independent and linear in the mass), 
while the light part scales with respect to topology, mass and volume essentially as 
predicted by chiral perturbation theory. The analysis in the chiral effective theory in 
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Figure 4: Results for {a^Xv) I Q-i^ as a function of the quark mass, for the lattice cl at 
several values of v. The data on the left are obtained with low-mode averaging, while 
on the right we show the values computed with the local estimator. 

the previous section suggests that these features are mostly volume independent. 

The effect of low-mode averaging can be appreciated in Fig. [U where we show 
{c?Xv)l0'm for v = 0-5 as a function of the quark mass, with (left) and without (right) 
LMAo. It is clear that the variance reduction is much more effective for the sectors with 
lower topological charge, which are dominated by infrared contributions. Nevertheless 
with LMA we still obtain a variance reduction of a factor ~ 2 up to = 5 in all our 
runs. 

5 Comparison with the effective theory 

In this section we compare our numerical results with the predictions of the chiral 
effective theory. We first focus on the finite-size scaling and the topology dependence of 
the condensate, and then we extract the value of Sefj from a fit in the mass of our data. 
Our raw numerical results of Xv for = 0-5, are reported in Table [2l The statistical 
errors have been estimated with a jackknife procedure. 

5.1 Finite-size scaling 

The parameters of the lattices cl and c2 have been chosen to carry out a finite-size 
scaling study. The lattice spacing is the same and the physical volume differ by more 
than a factor 3. The bare combinations a?{xi — X2)/clt^ and a?{x2 — Xs)/'^^'^ computed 

^At lower topologies and without LMA the average values and their errors are only indicative since 
the statistical analysis is invalidated by extreme statistical fluctuations. 
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lat 


u\am 


0.001 


0.003 


0.008 


0.012 


0.016 


cl 









0.00484(18) 


0.00610(16) 


0.00721(15) 




1 


0.000382(12) 


0.00114(3) 


0.00291(7) 


0.00422(8) 


0.00546(8) 




2 


0.000304(3) 


0.000910(8) 


0.002413(19) 


0.00359(3) 


0.00475(3) 




3 


0.000283(3) 


0.000849(8) 


0.002261(21) 


0.00338(3) 


0.00450(4) 




4 


0.0002712(22) 


0.000814(7) 


0.002168(17) 


0.00325(3) 


0.00433(3) 




5 


0.0002622(23) 


0.000787(7) 


0.002097(18) 


0.00314(3) 


0.00419(4) 


lat 


v\am 


0.000316 


0.000949 


0.00253 


0.00380 


0.00506 


c2 











0.0043(3) 


0.00480(25) 




1 






0.00182(12) 


0.00256(14) 


0.00319(15) 




2 


0.000157(9) 


0.00047(3) 


0.00123(6) 


0.00182(9) 


0.00237(10) 




3 


0.000129(3) 


0.000387(8) 


0.001026(20) 


0.00153(3) 


0.00202(4) 




4 


0.000110(4) 


0.000329(11) 


0.00088(3) 


0.00131(4) 


0.00174(5) 




5 


0.0001047(21) 


0.000314(6) 


0.000837(16) 0.001255(24) 


0.00167(3) 


lat 


u\am 


0.000612 


0.00184 


0.00490 


0.00735 


0.00980 


c3 









0.00246(11) 


0.00324(11) 


0.000394(10) 




1 


0.000217(7) 


0.000648(19) 


0.00168(4) 


0.00246(5) 


0.00321(5) 




2 


0.0001850(14) 


0.000556(4) 


0.001474(11) 


0.002200(15) 


0.002915(19) 




3 


0.0001729(18) 


0.000520(5) 


0.001382(14) 


0.002070(20) 


0.00276(3) 




4 


0.0001677(17) 


0.000504(5) 


0.001342(13) 


0.002012(20) 


0.00268(3) 




5 


0.0001622(17) 


0.000488(5) 


0.001299(13) 


0.001947(20) 


0.00260(3) 



Table 2: Numerical results for a'^Xv with v = 0-5. 



on these lattices are shown in the plots on the left of Fig. [5l Their values differ from zero 
by many standard deviations for both volumes. Finite size effects are clearly visible. 
The rescaled combinations are shown in the plots on the right of the same Figure as 
a function of {'mV)a~^ . The corresponding data sets are in very good agreement for 
all masses simulated within the statistical errors. These results are compatible with 
the finite size scaling behaviour expected at asymptotically large volumes for a theory 
which exhibits spontaneous symmetry breaking. Within our statistical errors we do not 
observe deviations from the leading scaling behaviour. Similar conclusions apply for the 
other combinations made with u = 0-3. Results from higher topological sectors tend to 
depart from the LO scaling behaviour, but our statistical errors are too large to draw 
any definite conclusion. 

5.2 First Leutwyler-Smilga sum rule 

The sensitivity on the topology of the spectral density of the low-lying eigenvalues in- 
duces a non-trivial i/-dependence in the chiral condensate. The latter has to be compared 
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Figures: The ratios a^(xi—X2)/(aW') and a^(x2— X3)/('^"^) (left), and a^(xi— X2)/("T'V') 
and a^(x2 — X3)/{^^) (right) as a function of (mV)a^^ for lattices cl and c2. 

with the prediction of the chiral effective theory in Eq. (j3.8p . For the ratios 



Xu3 XvA 

the prediction in the chiral limit is given in Eq. (j3.6p . It is parameter-free and valid up to 
the NLO. For the lattice cl and v = 1-5, the ratios in Eq. (|5.ip are shown in the first plot 
of Fig. [6] as a function of {mV)a~^ . As expected from Eq. (13. 8p the mass dependence 
is very mild, and the two lightest points are consistent with a flat behaviour for all 
combinations. Results for lattices c2 and c3 show analogous features. In the plot on the 
right we report the points at the lightest mass available for each of the three lattices, 
together with the theoretical expectations from the first Leutwyler-Smilga sum rule in 
Eq. (|3.6p . The very good agreement supports the fact that the topology dependence of 
these ratios is well reproduced in the (quenched) chiral effective theory. This is one of 
the main results of this paper. 

5.3 Discretization effects 

Since we use fermions with an exact chiral symmetry, the leading discretization effects 
in our observables are of order a?. We can estimate their magnitude by comparing 
renormalized dimensionless quantities calculated on the lattices cl and c3. In Fig. [7] we 
show the quantity mrQ{xui — Xi'2) ^ a function of mV/{ZsrQ). Within our statistical 
errors we do not observe significant deviations between the points of the two data 
sets. A similar behaviour is observed in the other topological sectors. Discretization 
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Figure 6: On the left, the ratios {xui — Xu2)/ixi ~ X'i) as a function of the quark mass 
for the lattice cl and for (1^1,1^2) = (1, 2), (1, 4), (1, 5). The points at m = represent 
the theoretical expectations from the Leutwyler-Smilga sum rule in Eq. ()3.6p . On the 
right the same ratios, at the lightest mass available for lattices cl, c2, c3, is compared 
with the theoretical expectations. 



effects are small and are likely to be comparable or below our statistical errors. This 
is not surprising: quenched computations of various physical quantities carried out 
with Neuberger fermions show small discretization effects at the lattice spacings of 
our simulations [13,54]. A continuum limit extrapolation of our results is beyond the 
scope of this paper, and we find it of limited interest before removing the quenched 
approximation . 

5.4 Extraction of the low-energy constant 

For every lattice and for {vi,i'2) = (0)1)) (1)2), (2,3), we extract the value of the 
low-energy constant T^^f^/Zs from a one-parameter fit of the observables Xui — Xu2i 
with the infrared tail truncated as described in Sec. HI to the corresponding functional 
form predicted by the chiral effective theory. Combinations with z^i , ^4 have larger 
statistical errors and are not considered in this analysis. The fit always reproduces the 
data very well, and the uncorrelated x^/d.o.f. are typically very small. The results for 
^es/Zs are reported in Tab. [3l At small values of (mV) the quenched effective theory 
predicts a logarithmic behaviour of the form (mV) log (mV) in XO) which is clearly seen 
in our data. If this term is removed from the fit function, the latter is not compatible 
with the data anymore. For illustration in Fig.[8]the function Xui ~Xu2 is superimposed 
on the data for Xui ~ Xu2 from lattice cl. 

The values of are in good agreement within each data set, a consequence 
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lat 




a^T,cs/Zs 


cl 


0-1 


0.0040(6) 




1-2 


0.0039(3) 




2-3 


0.0034(3) 


c2 


0-1 


0.0035(8) 




1-2 


0.0049(9) 




2-3 


0.0040(5) 


c3 


0-1 


0.0015(3) 




1-2 


0.00178(18) 




2-3 


0.00188(12) 



Table 3: Results for a'^Yl^g/Zs as determined from a fit of the lattice data (see text). 



of the fact that the effective theory reproduces the topology dependence observed in the 
data. The results from lattices cl and c2 confirm that, within our statistical errors, we do 
not observe a volume dependence in Sefj . By taking a = 0, the topological susceptibility 
from Ref. [57] and the quenched value of F from Ref. [29], the NLO formula in Eq. (I3.1ip 
suggests for Seff(i — 2.0 fm)/Scfr(-Z^ — 1.5 fm) a positive deviations from 1 of the order 
of 10 per cent. This is comparable to our statistical uncertainty, and the same holds if 
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Figure 8: The function Xui — Xu2 is superimposed on data for Xui — Xu2 from lattice cl. 



we vary a within a reasonable range of values Q Our best result for the renormalization 
group invariant condensate from the lattices cl and c3 is for the combination (i/i, V2) = 
(1,2), while from lattice c2 is for (2,3): 

(Seffrg)(L = 1.5 fm) = 0.33(3) cl , (5.2) 
(Seffrg)(L = 2.0 fm) = 0.34(5) c2 , (5.3) 
(Eeffrg)(L = 1.5 fm) = 0.29(3) c3 . (5.4) 

We take as our best estimate of the condensate at L = 1.5 fm the result from lattice 
c3, which is the one with the finer lattice spacing. The latter, converted into the more 
usual MS-scheme at 2 GeV [58], is 

r|jS^(2 GeV) = 0.40(4) at L = 1.5 fm, (5.5) 

and if we use the phenomenological value rg = 0.5 fm we obtain 

S^(2 GeV) = (290 ± 11 MeV)^ at L = 1.5 fm . (5.6) 

The errors in Eqs. (15. 5p and ()5.6p do not include uncertainties due to discretization 
effects. Our result is in the range expected from previous computations in the quenched 
approximation of QCD [5-14]. The determinations from the infinite volume regime [5,6, 
8,9] are affected by NLO chiral corrections which usually are not taken into account, and 
which are different from those in the finite volume regime. Even tough the comparison 

''Explicitly, one finds Ecft (L = 2.0 fm)/E<3ft (L = 1.5 fm) ~ 1 + 0.09 + 0.03a. 
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cannot be very accurate, it is reassuring the good agreement of our result with some 
of these determinations. First exploratory studies of the chiral condensate in the finite 
volume regime were performed with small volumes and low statistics [7, 10, 11]. This 
may be the reason why the large statistical fluctuations generated by the very small 
eigenvalues of the Dirac operator, which were first reported in Ref. [11], did not occur 
in some of these computations. A yet different method to extract was pursued in 
Ref. [12-14]. They compared the average value of the individual low-lying eigenvalues of 
the Dirac operator with the prediction of random matrix theory. Even if it relies on the 
non-trivial assumptions that the single eigenvalue of the Dirac operator renormalizes 
with XjZs and that the random matrix theory reproduces the QCD spectrum in the 
infrared, the values of Scfr obtained in these computations are in good agreement with 
our result. 

6 Conclusions 

The functional form of the (quenched) QCD chiral condensate in the volume range 
5-16 fm'^ and for masses mSF < 1 turns out to be in good agreement, within our 
statistical errors, with the prediction of the chiral effective theory: the volume, mass 
and topology dependence are well reproduced by our data. Close to the chiral limit we 
observe a detailed agreement with the first Leutwyler-Smilga sum rule. The low-energy 
constant Sefj is then extracted from a matching of the lattice results with the chiral 
formulas. 

A technical progress which made the computation possible is the introduction of 
a numerical estimator which is stable in the finite volume regime of (quenched) QCD. 
Very small eigenvalues of the Dirac operator, which do occur in the Monte Carlo history, 
do not generate large fluctuations, and a finite size scaling study of the condensate 
becomes feasible all the way down to the chiral limit. We have not tried to exploit the 
great freedom in the possible choices of the QCD observable, discussed in Sec. [21 which 
are unambiguously defined in the continuum limit and could lead to a more precise 
determination of the low-energy constant Seff . This is beyond the scope of this paper 
and we leave it to future studies. 

This technique is directly applicable to the computation of the condensate in full 
QCD simulations. 
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A Definitions and conventions 



We consider the Euclidean space-time discretized on a lattice with spacing a. Our 
conventions for the Dirac matrices are 

iji = In {in^ 71/ } = 2(5^1/ , (A.l) 

with 

75 = 70717273 = - 1 ) ■ ^^''^^ 

and the chiral projectors defined as P± = (1 it 75)/2. The massless Neuberger-Dirac 
operator is defined as 

^ = ^{l + 75sign(Q)} , (A.3) 

where 

Q = j5{aDw - 1 - s) , |s| < 1 , a = , (A.4) 

and Dw is the Wilson Dirac operator 

Dw = \ {7m(V; + V^) - aV;v4 . (A.5) 
The covariant forward and backward derivatives V;^ and V* are 

V^^(x) = ^ {Uix, fi)ilj{x + afi) - 4^{x)} , (A.6) 
V*^iPix) = ^{iPix)-U{x-afi,fi)-^iP{x-afl)}, (A.7) 

where U{x,fi) G SU(3) are the gauge variables and /t is the versor in direction /x. The 
massive quark propagator is 



2 

where 



Smix, y) = ^(x)V >(y) = ^ _ ^ jL*^ - - \ {x,y) , (A. 



= ( 1 - |Z) ) V . (A.9) 
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B Renormalization of the scalar density 



The bare singlet scalar density 

S^{x) =^{x)tP{x) (B.l) 

needs to be renormalized to make finite its correlation functions with other gauge in- 
variant operators inserted at a physical distance. itself and the identity multiplied 
by proper combinations of the mass matrix are the only operators with dimension d <3 
and the same transformations properties under the SU{Nf)i, (gi 5?7(A^/)r group. For 
quark masses degenerate and real, the renormalized operator can then be written as 

= Zs [S°{x) + cim + cam^l , (B.2) 

where Zs is the logarithmic-divergent renormalization constant of the scalar density, 
and ci and C2 are subtraction coefficients which at asymptotically large a diverge as 
1/a^ and ln(a) respectiveljU. The coefficients ci and C2 can be fixed, for example, by 
requiring that in the infinite volume limit the condensate satisfies 

{S') - = , (B.3) 

m=U 



dm 
dm? 



(5°) _ =0. (B.4) 



m=0 



The value of the renormalized condensate for m > is therefore prescription dependent, 
a well known fact in QCD. The cumulants of the topological charge Q are ultraviolet 
finite, and thus the distribution of the topological charge as defined in Eq. (12. 4p [18- 
20]. This implies that correlation functions of renormalized local operators inserted at 
a physical distance are finite also in the theory at fixed topology. In particular the 
combinations 

{S').. - {S\, (B.5) 

are unambiguously defined outside the chiral limit, i.e. they are independent on the 
particular prescription chosen to renormalize the chiral condensate at finite mass. 
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